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Abstract. Using ideal relativistic hydrodynamics in 2+1 dimensions, we study 
the collision energy dependence of radial and elliptic flow, of the emitted hadron 
spectra, and of the transverse momentum dependence of several hadronic parti- 
cle ratios, covering the range from Alternating Gradient Synchrotron (AGS) to 
Large Hadron Collider (LHC) energies. These calculations establish an ideal fluid 
dynamic baseline that can be used to assess non-equilibrium features manifest in 
future LHC heavy-ion experiments. Contrary to earlier suggestions we find that a 
saturation and even decrease of the differential elliptic flow V2(pt) with increas- 
ing collision energy cannot be unambiguously associated with the QCD phase 
transition. 

1 Introduction 

Relativistic hydrodynamics has proven to be very successful in describing the evolution of the 
hot and dense bulk matter created in heavy-ion collisions at the Relativistic Heavy Ion Collider 
(RHIC) |1I2I3I4| . Since, in the ideal fluid limit, it provides a direct connection between the 
equation of state (EOS) of the hot matter and the observed collective flow pattern that can 
be extracted from the measured hadron momentum spectra, this raised the hope of being able 
to constrain the QCD equation of state experimentally and to identify experimental signatures 
of the quark-hadron phase transition [5|6|7|8|9] which QCD predicts at high temperatures [TO] . 
Specifically, ideal fluid dynamics predicts a non-monotonic collision energy dependence of the 
PT-integrated elliptic flow t>2 [7], with a dip between SPS and RHIC energies caused by the 
softening of the QCD EOS around the quark-hadron phase transition. Unfortunately, large 
viscous effects during the late hadronic stage of the collision fireball expansion were found to 
spoil this phase transition signature |ll|12j , and experimentally the p^-integrated elliptic flow 
was found to increase monotonically with collision energy |13|14j (see also |15|16j ). 

Contrary to the p^-integratcd elliptic flow, however, the PHENIX Collaboration found that 
the differential elliptic flow V2(pt) for charged hadrons, at two fixed values of pr (j>t = 0.65 and 
1.75 GeV/c), does not grow monotonically with increasing collision energy, but instead appears 
to saturate at RHIC in the center-of-mass energy range between 63 and 200 GeV/nucleon pair 
|17j . In Refs. [18|19I20] this observation was brought into connection with the QCD quark- 
hadron phase transition, speculating that the non-monotonicity of elliptic flow caused by the 
softening EOS near the phase transition, in spite of being washed out in the p^-integrated V2, 
might survive viscous effects if measured at fixed pr- This speculation was based on the obser- 
vation that the monotonically increasing radial flow, which causes the transverse momentum 
spectra to fall off more slowly at higher collision energies, gives bigger weight to the high-p^ 
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region (where V2{pt) is larger than at low px) at high than at low collision energies, thereby 
counteracting in the /^-integrated elliptic flow any reduction of Vi^Pt) by a softening EOS. 
Any phase transition signature in the integrated v-i should thus manifest itself even more promi- 
nently in the differential elliptic flow vi(j>t) at fixed px, and it might thus remain visible in the 
energy dependence of V2(pt) even in the presence of viscous effects that smear out the signature 
in the pT-integrated v 2 |18I19| . 

This conjecture remained speculative as long as there existed no systematic calculation of 
the collision energy dependence of the hydrodynamically generated differential elliptic flow. 
The present work fills this hole in the literature. Our study is based on ideal fluid dynamics 
even though the ideal fluid assumption is known to gradually break down below RHIC energies 
|18j . We do not attempt to quantitatively describe relativistic heavy ion data, but to provide 
theoretical insights into the systematics of the energy dependence of radial and elliptic flow 
within the ideal fluid picture. However, since the validity of the ideal fluid picture is expected 
to improve with increasing collision energy, providing an almost quantitative description of 
most low-pT phenomena already at upper RHIC energies, the calculations presented here can 
be taken as a prediction for hadron spectra and their elliptic flow at LHC energies where effects 
from late hadronic viscosity are expected to become negligible [21]. As long as the viscosity 
of the quark-gluon plasma (QGP) stage (which at the LHC dominates over the hadronic one) 
remains sufficiently small, these predictions should give an accurate description of LHC datsQ- 

In addition to providing ideal fluid dynamical benchmarks for the LHC, an important finding 
of the present work is that a saturation and even decrease with growing beam energy of the 
differential elliptic flow V2{pr) at fixed pr cannot be unambiguously associated with the quark- 
hadron phase transition. We find that, at sufficiently high collision energies, V2(pr) at fixed 
Pt decreases with increasing beam energy even when the matter is initially so dense that 
all elliptic flow is generated far above the phase transition and therefore not affected by the 
softening of the EOS near the critical temperature T cr . This decrease of the pr-diffcrential 
elliptic flow at fixed pr is accompanied by a simultaneous increase of the pr-integrated elliptic 
flow. It is caused by an increase of the radial flow with growing beam energy which pushes 
the hadrons to larger pt and renders the momentum spectra less anisotropic at low pt- Our 
finding contradicts earlier speculations [L8I19I20] and makes the task of identifying the QCD 
phase transition experimentally even more difficult than previously thought. 

The present study is similar in spirit to and complements recent work by Nicmi and Eskola 
et al. who also made ideal fluid dynamical predictions for hadron spectra [35] and elliptic flow 
[23] at LHC energies. Our analysis goes beyond theirs in its systematic investigation of the 
beam energy dependence of these hadronic observables. 

2 Procedure 

Our hydrodynamic simulations were performed with the (2+l)-dimensional hydrodynamic code 
AZHYDRO that solves the equations of motion in the two directions transverse to the beam 
direction for an ideal fluid undergoing boost-invariant longitudinal expansion [7124125] . To run 
the simulations we must initialize the program by providing the following parameters: 

1. impact parameter b; 

2. initial proper time tq at which the fluid is considered to be in local thermal equilibrium and 
the hydrodynamic expansion stage begins; 

3. initial peak entropy density sq in central (b = 0) collisions: this parameter is used to control 
the final charged hadron multiplicity dN^/dy and serves as a proxy for the collision energy 
(see discussion below); 

1 As the beam energy and thus the initial fireball temperature T increases, the effective coupling 
strength u B (T) decreases logarithmically; hence, the shear viscosity to entropy ratio n/s is expected to 
slightly increase from RHIC to LHC, suggesting stronger dissipative effects at the LHC. However, the 
particle density n ~ T s and mean free path A ~ 1/T vary much more rapidly with T; correspondingly, 
the sound attenuation length _T S = n/(sT) decreases from RHIC to LHC, and at any given time r 
viscous effects, characterized by the ratio = f tj^, should be weaker at LHC energies than at RHIC. 
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4. peak value ub,o for the initial net baryon number density in central collisions: this parameter 
is used to control the finally observed anti-proton/proton ratio; 

5. equation of state (EOS): we use an EOS that matches an ideal gas of quarks and gluons at a 
critical temperature T cr = 164 MeV to a realistic hadron resonance gas with non-equilibrium 
hadron abundances whose chemical composition is frozen in chemical equilibrium at T cr 
[24126127] : 

6. freeze-out temperature T^ ec : this parameter controls where the hydrodynamic stage ends 
and how much collective flow builds up during the expansion. We use Td oc = 100 MeV. 

We simulate Au+Au collisions. The initial transverse entropy and baryon number density pro- 
files are taken to be proportional to each other (constant entropy per baryon) and are calculated 
from the Glauber model with a mixture of soft and hard collisions, assuming 75% weight for the 
soft contribution, proportional to the transverse density of wounded nucleons, and 25% weight 
for the hard contribution, proportional to the transverse density of binary nucleon-nucleon 
collisions (see [5] for details). This mixture produces the correct centrality dependence of the 
charged hadron multiplicity at RHIC [2], and we assume that it doesn't change appreciably 
between RHIC and LHC. For the systematics of our study this is not a critical assumption. 
The initial transverse flow velocity at To is assumed to vanish. 

In the absence of shocks, ideal fluid dynamics conserves entropy. Therefore, the observed 
charged hadron multiplicity dN^/dy, which measures the entropy of the final state, can be 
directly related to integral over the initial entropy density profile, parametrized by sq. Hy- 
drodynamics cannot predict its own initial conditions, and a fundamental theory that reliably 
predicts the charged multiplicity as a function of collision energy y/s does not exist yet. We 
therefore use dN^/dy or, equivalcntly, the initial peak entropy density in central Au+Au col- 
lisions, sq, as a proxy for the collision energy: At any given collision energy, a measurement 
of dN c h/dy in the most central collision events fixes the value of sq to be used in ideal fluid 
simulations at that energy. 

Using linear longitudinal expansion without transverse flow at very early times, dN^/dy 
and sq are thus related by 



dN ch 
dy 



t J d 2 x± s(x±,t q ) oc t s - 



cx t / d 2 x ± s(x±,t q ) oc t s - (1) 



This relation involves two constants of proportionality, of which the first is the inverse of the 
average entropy per charged particle in the final state and is thus determined by the EOS at 
freeze-out, while the second depends on the shape of the initial entropy density profile which 
is fixed by the employed Glauber model. Equation (TT]) shows that s and t are inversely 
proportional to each other. Since we vary the collision energy over orders of magnitude, the 
initial energy and entropy densities likewise vary by large factors. It is therefore not adequate 
to assume a constant thermalization time To independent of collision energy. Due to larger 
densities and temperatures at higher collision energies, thermalization should happen more 
rapidly. We assume that To scales with collision energy by satisfying a quantum mechanical 
uncertainty relation between the initial temperature To (i.e. the initial average thermal energy) 
and the initial time To: To To = const. [28j . The constant is fixed by assuming the standard value 
To = 0.6fm/c at RHIC energies, y 7 ! = 200 GcV/nucleon pair (it turns out to be approximately 
1 [2]). Since for all collision energies considered here the fireball center in b = collisions is in 
the QGP phase, sq cx Tq . Combining these relations with Eq. (fT]), one finds 



tq (y/si) _ T (ygj) _ I dN ch / dy (y^) s (ygl) _ / dN ch /dy (y/s!) \ 1 

to (y/sj) ~ T (VST) ~ V dN ch /dy (^l) ' so (yW) ~ \dN ch /dy (y^) J 



(2) 



Figure[T]shows a compilation [25] of measured hadron multiplicities from (A«200) + (Aw200) 
collisions at a variety of collision energies explored at the AGS, SPS and RHIC, together with 
a linear fit in log 10 y^: 

^ = 312.5 log 10 V^- 64.8 (3) 
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Fig. 1. (Color online) Charged particle multiplicity dN^/dr] vs. center of mass energy per nucleon pair 
y/s (lower horizontal axis) and initial peak entropy density so (upper horizontal axis). The experimental 
data are from central Au+Au and Pb+Pb collisions at AGS, SPS, and RHIC energies and were compiled 
in Ref. [29]. The dashed line is a linear fit to the data. Crosses indicate values for so and dN c h/dr] for 
which hydrodynamic simulations were performed. See text for details. 



The upper horizontal axis maps the multiplicities on the vertical axis onto sq values, using 
Eq. © with s (\/i=200AGeV) = 117/fm 3 [25l2] , The linear fit © provides guidance for 
which values to expect for dN^/drj and s a t hitherto unexplored LHC energies. We emphasize, 
however, that So is not directly related to y/s, but only indirectly via Eq. ^ through the charged 
multiplicity measured at that value of y/s. For this reason we have performed several "LHC 
simulations" with a variety of s values within a reasonable range of the value predicted by the 
simple-minded linear extrapolation of existing data shown in Fig. [1] The first day of LHC heavy- 
ion experiments will tell us which of these simulations to choose as hydrodynamic reference for 
the data. 



We finally comment on our choice of the initial peak net baryon density nsfi- At RHIC 
energies it was fitted to the measured p/p ratio to give i%b,o = 0.44/fm 3 at To = 0.6fm/c. We 
did not attempt to fit the -y/s-dependence of ns,o from p/p ratios measured at lower collision 
energies. Instead we simply held Hb,o = 0.44/fm 3 constant over the entire energy range while 
reducing To with increasing collision energy according to Eq. This implies that, at fixed 
r, tib decreases with increasing multiplicity as 1/ y/dN^/dy. Together with the increasing sq 
values, this leads to a significant (although perhaps not quite strong enough) decrease of the 
net baryon to entropy ratio at midrapidity. Correspondingly, the baryon chemical potential at 
chemical decoupling decreases, reflecting increasing nuclear transparency and decreasing baryon 
stopping power at higher energies. At our highest sq value, so = 271/fm 3 (last cross in Fig.[T]), 
we simply set ns o — 0, assuming approximate baryon-antibaryon symmetry at midrapidity at 
the LHC. 
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Fig. 2. (Color online) Transverse momentum spectra for thermal pions (n + , upper panels) and protons 
(p, lower panels) at low (left panels) and intermediate (right panels) values of pr, for central Au+Au 
collisions with a variety of initial peak entropy densities so- 



3 Radial Flow 



3.1 Transverse momentum spectra from AGS to LHC 

Figure [2] shows the p^-spectra for directly emitted ir + mesons (upper row) and protons (lower 
row). In this plot we have neglected feed-down from resonance decays; its inclusion is com- 
putationally intensive but, since we keep T^ ec the same at all collision energies, it will not 
qualitatively affect the systematics shown in Fig. [21 Here and in the following the curves are 
labelled by the value of sq. The reader can use Fig. Q] to translate this value into charged hadron 
multiplicities (which include resonance feeddown) and to estimate the corresponding collision 
energy. 

The normalization of the pion p^-spectra is seen to increase with so, reflecting the growth of 
the total multiplicity with increasing collision energy. They also become systematically flatter 
as so increases; since the decoupling temperature Td oc is held fixed, this is an unambiguous 
signature for larger radial flow at higher collision energies. 

The radial flow effects on pion and proton spectra are similar at large transverse momenta 
(right pnaels in Fig. [2]) where their rest mass difference can be neglected and the flow effect 

can be understood in terms of the simple blucshift formula T s io P o = ^decy^j^j^} [30131] . At 
low pt, however, they are much more pronounced for the heavier protons (lower left panel in 
Fig. [2]) |30I31| . One sees that the radial flow pushes the protons away from pt = 0, leading to 
both a dramatic flattening of the pT-spcctrum and a decrease of the proton yield at low pt, in 
spite of the overall increase in multiplicity. (It should be mentioned that the decrease in proton 
yield at low pt is amplified by the decrease of the baryon chemical potential at higher collision 
energies, but we checked that it is also visible in the antiproton spectra.) Only at the highest 
collision energies, where the low-pr proton spectra are almost flat, the low-p^ proton yields are 
seen to follow the general increase in multiplicity. 
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3.2 px- and m^-dependence of particle ratios: RHIC vs. LHC 



Hydrodynamic radial flow, which leads to flatter p^-spectra for heavy particles, is a key con- 
tributor to the observed [4] strong rise of the p/ir and A/K ratios at low pr at RHIC [2]- The 
left column in the left panel of Figure [3] shows that this rise is predicted to be slower at the 
LHC than at RHIC. Since all spectra are flatter at the LHC due to increased radial flow (right 
column of the left panel) while their asymptotic ratios at pt — * oo (given by their fugacity and 
spin degeneracy ratios [5]) remain similar, the ratios between the spectra increase more slowly 
with pt- 




p T (GeV) p T (GeV) m T -m (GeV) ro^GeV) 



Fig. 3. (Color online) Left panel: Transverse momentum spectra (right half) and selected particle 
ratios (left half) as functions of pr for central Au+Au collisions at RHIC and LHC energies (see text 
for details). A zoomed version of this plot which focusses on the region pr < 2GeV/c can be found in 
Fig. 53 of Ref. [32] , Right panel: Same as left panel, but plotted as a function of transverse mass mr or 
transverse kinetic energy mr — mo, respectively, over the mr range where hydrodynamics is expected 
to be a valid description. All curves shown in this Figure include the contributions from resonance 
decays. 



In the right panel of Fig. [3] we redraw the curves shown in the left panel as a function of 
transverse mass mj- (for the spectra) or transverse kinetic energy m^-mo (for the particle 
ratios). We do so in order to isolate flow effects, by eliminating the kinematic contribution 
to the rise of the pT-dependent heavy/light hadron ratios that results from plotting thermal 
distributions (which depend on jr 1 -) as a function of px = y m^p—m^. In the absence of 
radial flow, the my-spectra would show perfect TOr-scaling, and (except for minor effects from 
resonance feeddown and Bose statistics for pions) the ratios between spectra of particles with 
different rest masses would thus be independent of tot or tot — too. The rise of the ratios shown 
in the left half of the right panel of Fig. [3] (which is much weaker than that in the left half of 
the left panel) can thus be attributed almost exclusively to radial flow effects. A comparison of 
future LHC data on (tot— TOo)-dcpendcnt heavy/light particle ratios with the hydrodynamical 
predictions shown in Fig. [3] can therefore help to separate collective flow effects from more 
exotic explanations such as "baryon junctions" [32] • 
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Fig. 4. (Color online) Differential elliptic flow V2(pr) for directly emitted pions (tt + , upper panels) and 
protons (p, lower panels) at low (left panels) and intermediate (right panels) values of pt, for b = 7fm 
Au+Au collisions with a variety of initial peak entropy densities so. 

4 Elliptic Flow 

While ideal fluid dynamics begins to break down below RHIC energies, due to viscous effects in 
the late hadronic stage which persist even at RHIC |11I12| , its validity is expected to improve 
at the LHC where the elliptic flow saturates in the QGP stage and effects from late hadronic 
viscosity become negligible [21]. Early viscous effects in the QGP stage seem small at RHIC 
[2|12j . and recent results from Lattice QCD indicate little change of its specific shear viscosity 
rj/s from RHIC to LHC [33]. The following ideal fluid dynamical predictions for soft (jpr < 
2-3GcV/c) hadron production in (^4«200)+(A«200) collisions at the LHC should thus be 
quite robust. 

In Figure [4] we look at the differential elliptic flow for both protons and ir + as a function 
of transverse momentum for b = 7fm Au+Au collisions at various collision energies. The left 
panels are expanded views of the low-px end. At all collision energies, over 99 percent of all 
particles are emitted with transverse momenta below 1.5 GeV/c; for this reason we focus our 
attention on the system's characteristics at low pt0 

By scanning the curves in order of increasing Sq one notices for both pions and protons 
that, at low pt, the differential elliptic flow at fixed pt is not monotonic with collision energy. 

2 The astute reader may be puzzled, as we were initially, by the fact that, for the highest collision 
energies studied here, the elliptic flow peaks at intermediate pt and then decreases again, instead 
of monotonically increasing with pr- This effect appears to be caused by the following phenomenon: 
even though on average the hydrodynamic flow at freeze-out is stronger in the reaction plane than 
perpendicular to it, at high collision energies the largest value of the flow velocity is found on the 
freeze-out surface in out-of-plane direction. Since the highest-pr hadrons are emitted from fluid cells 
with the largest flow velocities, this causes a reduction of vi at very high px, causing it eventually to 
even turn negative [23]. We thank P. Huovinen for a clarifying discussion on this point. 
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Fig. 5. (Color online) Upper panel: pr-integrated elliptic flow for pions from b = 7fm Au+Au collisions 
as a function of collision energy, assuming the relation shown in Fig. [TJbetween so and i/s. The solid line 
takes only directly emitted thermal pions into account, the dashed line includes pions from resonance 
decays. Lower panel: Energy dependence of thermal pion elliptic flow from b = 7 fm Au+Au collisions 
at four fixed values of pr (in GeV/c), as indicated. 



At the low-energy end V2(pt) first increases with so, but when the collision energy is further 
increased beyond yfs > 10 GeV/nucleon pair, the differential elliptic flow decreases again. 

This is explored more quantitatively in the lower panel of Fig. \5\ which shows the differential 
elliptic flow of directly emitted pions at fixed pr as a function of ^fs (using the fit in Fig. Q] 
to translate Sq into \/s), f° r f° ur different pr values. The three curves corresponding to pr = 
0.1, 0.5 and 1 GcV/c all show this non-monotonic behavior of a V2(pt) that first rises and then 
decreases with collision energy. The increase at low y/s values is easy to understand: at low 
collision energies, the fireball decouples before the elliptic flow can saturate; by increasing the 
collision energy, the fireball lifetime is increased, allowing the elliptic flow to grow and reach 
values closer to its asymptotic saturation value. Between SPS and RHIC energies, the fireball 
lifetime is sufficiently long for the elliptic flow to more or less saturate [7]. The explanation for 
the decrease of V2(pt) at higher energies is more subtle: As the collision energy increases further, 
the magnitude of the radial flow continues to grow, too, leading to ever flatter p-r-spectra. 
Flatter p^-spectra show less azimuthal variation and thus exhibit smaller V2 coefficients [35] . 
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Fig. 6. (Color online) Differential elliptic flow V2(pt) for thermal pions and protons from b = 7fm 
Au+Au collisions at RHIC (s = 117 fin" 3 ) and LHC (s = 271 mr" 3 ) energies, plotted together for 
comparison. 

This flattening effect of radial flow on the j>T-spectra is particularly pronounced at low px (see 
Fig. [2]). This may explain why the non-monotonic beam energy dependence of V2(px) is only 
seen at low px up to about 1 GeV/c and seems to disappear above px = 2 GeV/c (dashed line 
in the lower panel of Fig. [5]). We expect the non-monotonicity of V2(pt) to persist to larger pt 
values for protons whose Iow-pt spectra are more strongly affected by radial flow than those of 
pions (see Fig.[2|). Figure [6] supports this expectation by showing that for protons (in contrast to 
pions) the differential elliptic flow at the LHC is smaller than at RHIC over the entire p*r-range 
shown in that Figure. 

It is important to note, however, that, while the differential elliptic flow at fixed px de- 
creases between RHIC and LHC energies, the total momentum anisotropy, reflected by the 
PT-integratcd elliptic flow coefficient V2, increases from RHIC to LHC. Radial flow flattens the 
LHC spectra dramatically, putting a larger weight on the larger V2 values at higher px, and 
as a result the px -integrated elliptic flow is larger at the LHC, even though at any fixed px 
value below 1 GeV/c it is smaller than at RHIC. This is shown in the upper panel of Fig. [5] for 
both directly emitted pions (dashed line) and all pions including resonance decay contributions 
(solid line). Again, one sees first a rise of V2 at low yfs (AGS energies), followed by a decrease 
between SPS and RHIC and another rise between RHIC and LHC [7]. The non-monotonicity 
of the px -integrated V2 can be related to the quark- hadron phase transition and the corre- 
sponding softening of the EOS in the transition region [7]. In the experimental data this phase 
transition signature is, unfortunately, washed out by strong viscous effects in the late hadronic 
stage of the fireball expansion which become increasingly important at low collision energies 
[11112118] (where the fireball spends most of its time in the hadronic phase). As a result, the 
experimentally measured integrated elliptic flow V2 rises monotonically with yfs, approaching 
the ideal fluid limit only at or above RHIC energies [21] . 

The upper panel of Fig. [5] shows that at all collision energies the directly thermally emitted 
pions show more p^-integrated elliptic flow than all pions together. From this we can deduce 
that pions emitted from resonance decays have a lower momentum anisotropy. This is at least 
partially understood by the fact that decay pions typically have smaller transverse momenta 
than their parent resonances [31) . and that the parent resonances, being heavier, contribute less 
elliptic flow at low px than pions (see Fig. [6]) [36] . 

Comparison of the curves shown in the upper panel of Fig. [5] with those in the right panel 
of Fig. 2 in Ref. [21] shows that in a hybrid model approach (which includes viscous effects 
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in the hadronic phase) the p-r-integrated elliptic flow increases by about 25% between RHIC 
and LHC whereas in the ideal fluid approach (which neglects hadronic viscosity) it increases 
by less than 10%. The largest contribution to the expected increase of from RHIC to LHC 
is therefore due to the disappearance of late hadronic viscous effects between RHIC and LHC 
because the fireball spends a smaller fraction of it time in the hadronic phase and the elliptic 
flow saturates before the QGP converts to hadrons. 



5 Conclusions 

We have shown that, while the radial flow and the pr-integrated elliptic flow increase from 
RHIC to LHC, the differential elliptic flow at fixed pr < 1.5GeV/c decreases in the same 
collision energy range. This decrease of t>2(pr) is driven by a (relative) depletion of low pt 
hadrons by radial flow which pushes the hadrons to larger transverse momenta. The observed 
"saturation" of ^(pt) seen by PHENIX [T7] between y/s = 63 and 200 GeV likely signals the 
onset of this kinematic effect and should thus not be interpreted without further scrutiny as a 
signature for the quark-hadron phase transition. 

When experimental data will become available for very high energy collisions at the LHC, 
the ideal fluid dynamical calculations presented here can serve as a benchmark for comparison 
with experiment, providing insights into the validity of ideal hydrodynamics at high energy 
densities and permitting quantification of deviations from ideal fluid behaviour. This should 
facilitate the discovery of possible novel effects, beyond those expected from collective dynamics, 
at LHC energies. 
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